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Abstract 

Introduction: The role of interleukin (IL)-lp -31T/C promoter polymorphism 
in the pathogenesis of chronic obstructive pulmonary disease (COPD) has 
been studied with inconsistent results. This meta-analysis was performed to 
assess the association of IL-ip -31T/C promoter polymorphism with COPD 
susceptibility. 

Material and methods: Published case-control studies from PubMed and 
China National Knowledge Infrastructure (CNKI) databases were retrieved. 
Data were extracted and pooled odds ratios (OR) with 95% confidence inter- 
vals (CI) were calculated. 

Results: Six case-control studies were included in this meta-analysis. The pooled 
effect size showed that I L- 1 p -31T/C was significantly associated with COPD 
susceptibility in an overdominant genetic model (CC+TT vs. TC, OR: 0.77, 
95% CI: 0.63-0.94), indicating that homozygotes (CC and TT) had a de- 
creased risk for COPD compared with heterozygotes (TC). In the subgroup 
analysis by ethnicity, the results indicated that I L- 1 p -31T/C was signifi- 
cantly correlated with COPD susceptibility in Asians (overdominant model, 
OR: 0.75, 95% CI: 0.61-0.93), further suggesting a protective role of IL-lp 
-31T/C in COPD pathogenesis in Asians. Moreover, after excluding the study 
without Hardy-Weinberg equilibrium, the pooled results were robust and no 
publication bias was found in this study. 

Conclusions: This meta-analysis suggests that IL-lp -31T/C promoter poly- 
morphism confers protection against COPD in Asians. 

Key words: chronic obstructive pulmonary disease, interleukin-ip, poly- 
morphism. 



Introduction 

Chronic obstructive pulmonary disease (COPD) is a worldwide dis- 
ease, characterized by persistent airflow limitation that is progressive 
and associated with an enhanced inflammatory response in airways and 
lungs [1]. To date, although the underlying mechanisms of COPD remain 
not fully elucidated, a genetic predisposition of COPD has been strongly 
evidenced [2-4]. 

Interleukin (IL)-1|3, a proinflammatory cytokine, is encoded by the 
I LI B gene on chromosome 2ql4, which is identified as a candidate gene 
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of COPD [5]. In the past decade, the relationship 
of polymorphisms in the IL1B gene with COPD 
risk has been well studied, including the -31T/C 
promoter polymorphism in the IL1B gene with 
the potential to alter the transcriptional activity 
[6]. However, inconsistent results were reported 
[7-11], which may be owing to differences in eth- 
nicity and sample size in individual studies, result- 
ing in lower statistical power. It is suggested that 
meta-analysis could be a useful means to pool the 
independent statistical powers and thus achieve 
a quantitative understanding of the associations. 
Accordingly, in the present study, a meta-analysis 
was performed to determine IL-1|3 -31T/C pro- 
moter polymorphism and the risk of COPD. 

Material and methods 

Search strategy 

The literature search was conducted using 
PubMed and China National Knowledge Infra- 
structure (CNKI) (http://www.cnki.net/) databas- 
es. CNKI database was founded by Tsinghua Uni- 
versity of China in 1996 and includes over 8000 
Chinese journals covering natural and social sci- 
ences. The languages were limited to English and 
Chinese. The following search terms were utilized: 
interleukin lbeta or interleukin lb or ILlbeta or 
ILlb, and polymorphism or variant or SNP or gen- 
otype, and chronic obstructive pulmonary disease 
or COPD. 

Data extraction 

Two independent reviewers collected the data 
according to the inclusion and exclusion criteria. 
For inclusion in the meta-analysis, retrieved arti- 
cles had to inform about the number of cases and 
controls, and the number of individual genotypes 
in cases and controls. Exclusion criteria in the me- 
ta-analysis were 1) not case-control genetic study, 
2) duplicated report, 3) no useful data reported, 
4) other IL-1(3 polymorphisms than -31T/C. Un- 
published data were not considered. Disagree- 
ment was resolved by discussion before reaching 
a consensus. If more than one article was pub- 
lished by the same group using the same cases, 
the study with a higher sample size was selected. 

Statistical analysis 

Categorical variables are presented as odds ra- 
tio (OR) with 95% confidence interval (CI). CC, TC 
and TT are the genotypes of I L- 1 (3 -31T/C polymor- 
phism. OR1, OR2 and OR3 were calculated as fol- 
lows: OR1: CC vs. TT; OR2: TC vs. TT; OR3: CC vs. 
TC. These pairwise differences (OR1, OR2 and OR3) 
were used to indicate the most appropriate genetic 
model as follows: if OR1 = OR2 * 1 and OR3 = 1, 
then a dominant model was suggested; if OR1 = 



OR3 ^ 1 and OR2 = 1, then a recessive model was 
suggested; if OR2 = 1/OR3 * 1 and OR1 = 1, then 
a complete overdominant model was suggested; if 
OR1 > OR2 > 1 and OR1 > OR3 > 1 (or OR1 < OR2 
< 1 and OR1 < OR3 < 1), then a codominant model 
was suggested. Once the best genetic model was 
identified, this model was used to collapse the 
three genotypes into two groups (except a codom- 
inant model) and to pool the results again. Pooled 
ORs with 95% CI were calculated and p < 0.05 was 
accepted as statistical significance. Heterogeneity 
was checked by the Q test. Meta-analysis was done 
with the fixed-effects model when there was no 
heterogeneity (p > 0.1). Otherwise, the random-ef- 
fects model was used. Subgroup analysis was per- 
formed by ethnicity to assess the effect of possible 
clinical heterogeneity on the summary ORs. Pear- 
son's x 2 test was use d to determine whether 
the observed frequencies of genotypes in controls 
conformed to the Hardy- Weinberg equilibrium 
(HWE). Studies with controls that depart from HWE 
(p < 0.05) were subjected to a sensitivity analysis in 
order to check the consistency of the overall effect 
size. Funnel plots, as well as Begg's rank correlation 
test and Egger's linear regression test, were used to 
assess the potential publication bias, and p < 0.05 
was considered significant publication bias. All 
analyses were conducted using Revman 5.0 (Ox- 
ford, UK, The Cochrane Collaboration) and Stata 
11.0 (StataCorp LP, College Station, TX, USA) [12]. 

Results 

Studies included in the meta-analysis 

Sixteen studies were relevant to the search 
terms. After reviewing the titles, abstracts and 
articles, 11 studies were excluded; thus 5 articles 
with 6 case-control studies for -31T/C matched 
the inclusion criteria (Figure 1). Of the 5 included 



16 articles relevant to the search terms 



9 articles excluded after title and abstract review: 
6 not case-control genetic studies 
3 other IL-ip polymorphism studies 



7 articles for full-text review 



2 artcles excluded after full-text review 
2 no useful data reported 



5 articles included in this meta-analysis 



Figure 1. Flow diagram of search process 
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Table I. Main characteristics of included studies 



Ref. 


Country 


Race 


Genotyping 


Source of 
control 




COPD 








Control 














Total 


CC 


TC 


TT 


Total 


CC 


TC 


TT 


HWE (p) 


~w~ 


Japan 


Asian 


PCR + 
RFLP 


Healthy 
population 


68 


20 


32 


16 


85 


20 


40 


25 


0.8771 


[5] 


Japan 


Asian 


PCR + 
RFLP 


Healthy 
population 


88 


16 


52 


20 


60 


8 


31 


21 


0.8082 


[5] 


Egypt 


Caucasian 


PCR + 
RFLP 


Healthy 
population 


105 


11 


45 


49 


71 


16 


29 


26 


0.3734 


[6] 


Korea 


Asian 


ABI 
sequencer 


Healthy 
men 


311 


74 


179 


58 


386 


109 


177 


100 


0.2696 


[7] 


China 


Asian 


PCR + 
RFLP 


Healthy 
population 


162 


29 


101 


32 


162 


27 


99 


36 


0.0159 


[8] 


Taiwan 


Asian 


PCR + 
RFLP 


Healthy 
population 


30 


6 


18 


6 


115 


23 


64 


28 


0.4665 



H WE - Hardy-Weinberg equilibrium, PCR - polymerase chain reaction, Ref- reference, RFLP - restriction fragment length polymorphism 



articles, 4 were published in English, and 1 was 
published in Chinese. These studies were carried 
out in China, Japan, Korea, Egypt, and Taiwan. No- 
tably, the study by Hegab et al. [8] was performed 
in two ethnic populations (Japan and Egypt). The 
main features of the studies included in this me- 
ta-analysis are presented in Table I. 

Quantitative synthesis 

In order to indicate the most appropriate ge- 
netic model, OR1, OR2 and OR3 were calculated. 



Results showed OR1 = 1.12, OR2 = 1.38 and 
OR3 = 0.80 for -31T/C, suggesting an overdom- 
inant genetic model (CC+TT vs. TC). The pooled 
effect size showed that there was a significant 
association of IL-ip -31T/C with the risk of COPD 
(OR: 0.77, 95% CI: 0.63-0.94, p = 0.009, fixed 
model, Figure 2). In the subgroup analysis by eth- 
nicity, the results indicated that I L- 1 p -31T/C was 
significantly correlated with COPD susceptibility in 
Asians (OR: 0.75, 95% CI: 0.61-0.93, p = 0.008, 
fixed model, Figure 3). 



Studies COPD (n//V) Control (n/W) Weight [%] OR (95% CI) 



Asada etal. 2005 


36 


68 


45 


85 


8.6 


1.00 (0.53, 1.89) 






Hegab etal. 2005 


36 


88 


29 


60 


9.3 


0.74 (0.38, 1.43) 






Hsieh etal. 2008 


12 


30 


51 


115 


5.8 


0.84 (0.37, 1.90) 






Lee et al. 2008 


132 


311 


209 


386 


48.8 


0.62 (0.46, 0.84) 






Liu etal. 2012 


61 


162 


63 


162 


17.9 


0.95 (0.61, 1.49) 






Hegab etal. 2005 


60 


105 


42 


71 


9.8 


0.92 (0.50, 1.70) 






Total (95% CI) 




764 




879 


100.0 


0.77 (0.63, 0.94) 


♦ 




Total events 


337 




439 












Heterogeneity: % 2 = 


3.72, df 


= 5(p = 


0.59), I 2 


= 0% 




1 


—\ 


i 1 1 



Test for overall effect: Z= 2.61 (p = 0.009) 0.02 0.1 1 10 50 

Test for subgroup differences: y 2 = 0.39, d/= 1 (p = 0.53), I 2 = 0% Decrease risk Increase risk 

Figure 2. Forest plots of OR with 95% CI for the association of IL-ip -31T/C and COPD risk (CC+TT vs. TC) 



Studies COPD (n//V) Control (n/W) Weight [%] OR (95% CI) 



Asada etal 2005 


36 


68 


45 


85 


8.6 


1.00 (0.53, 1.89) 






Hegab etal. 2005 


36 


88 


29 


60 


9.3 


0.74 (0.38, 1.43) 






Hsieh etal. 2008 


12 


30 


51 


115 


5.8 


0.84 (0.37, 1.90) 






Lee et al. 2008 


132 


311 


209 


386 


48.8 


0.62 (0.46, 0.84) 






Liu etal. 2012 


61 


162 


63 


162 


17.9 


0.95 (0.61, 1.49) 






Total (95% CI) 




659 




808 


90.2 


0.75 (0.61, 0.93) 


♦ 




Total events 


277 




397 






1 


— 1 


1 1 



Heterogeneity: y 2 = 3.33, df= 4 (p = 0.50), I 2 = 0% 0.02 0.1 1 10 50 

Test for overall effect: Z = 2.67 (p = 0.008) Decrease risk Increase risk 



Figure 3. Forest plots of OR with 95% CI for the association of IL-lp -31T/C and COPD risk in Asians (CC+TT vs. TC) 
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Test of heterogeneity 

Significant heterogeneity was not revealed be- 
tween all studies in CC+TTvs. TC comparison (over- 
dominant model) for -31T/C {I 2 = 0%, p = 0.59). 
When stratified by ethnicity, no heterogeneity was 
observed between the studies in Asians. 

Sensitivity analyses 

In the present meta-analysis, in only one study 
[10] was there lack of HWE, which had a potential 
to influence the robustness of the present me- 
ta-analysis. After exclusion of this study, the pat- 
tern of the pooled effect size persisted, while the 
overall heterogeneity was not significantly altered 
for -31T/C {I 2 = 0%, p = 0.62), indicating that this 
study might not contribute to the overall hetero- 
geneity. 

Publication bias 

The funnel plots showed no significant asym- 
metry in this meta-analysis (Figure 4). Moreover, 
publication bias was not suggested by Begg's 
rank correlation test [p = 1.000) or Egger's linear 
regression test (p = 0.138). 

Discussion 

IL-ip, a proinflammatory cytokine, contributes 
to the pathogenesis of COPD. Recently, associa- 
tion of IL-1(3 -31T/C promoter polymorphism with 
COPD risk caught more attention, owing to its po- 
tential to alter the transcriptional activity. 

The pooled effect size showed that IL-ip 
-31T/C was significantly associated with COPD 
susceptibility in an overdominant genetic model, 
indicating that homozygotes (CC and TT) had a de- 
creased risk for COPD compared with heterozy- 
gotes (TC). Significant heterogeneity was not re- 
vealed between all studies in CC+TT vs. TC com- 
parison for -31T/C. In the subgroup analysis by 
ethnicity, the results indicated that IL-ip -31T/C 
was significantly correlated with COPD suscep- 
tibility in Asians, further suggesting a protective 
role of IL-ip -31T/C in COPD pathogenesis in 
Asians. After exclusion of the study without HWE 
by Liu et al. [10], the pattern of the pooled effect 
size persisted, while the overall heterogeneity was 
not significantly altered for-31T/C, indicating that 
this study might not contribute to the overall het- 
erogeneity. Publication bias was not suggested in 
the present study, possibly owing to the deliberate 
search strategy and data extraction. 

However, some limitations should be consid- 
ered in this meta-analysis. First, large sample size 
studies were lacking and the pooled statistical 
power might be insufficient. Second, the pooled 
estimates were not based on adjustment by con- 




-l - 



0 0.2 0.4 

s.e. of logor 

Figure 4. Begg's funnel plot for evaluation of publi- 
cation bias in the included studies on the associa- 
tions of IL-ip -31T/C with COPD risk 

founding factors, such as sex, age, and smoking 
history. Third, the lack of original data in the stud- 
ies limited our further analysis of the potential in- 
teractions between gene and gene, or gene and 
environment, which might modulate COPD risk. 

In conclusion, although the pooled estimates 
should be interpreted with caution, our meta- 
analysis suggests that IL-ip -31 T/C promoter 
polymorphism confers protection against COPD 
in Asians. However, large sample size studies with 
unbiased genotyping methods, standardized de- 
fined COPD cases and matched controls in dif- 
ferent populations, as well as more detailed data 
about individual and environment, are warranted. 
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